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O 

CN ! Abstract 

> 
O 

I We show that models for screening of nuclear reactions in the Sun can be 

\^ ■ tested by means of helioseismology. As well known, solar models using the 

weak screening factors are in agreement with data. We find that the solar 
. model calculated with the anti screening factors of Tsytovitch is not consis- 

\ tent with helioseismology, both for the sound speed profile and for the depth 

of the convective envelope. Moreover, the difference between the no-screening 
and weak screening model is significant in comparison with helioseismic un- 
certainty. In other words, the existence of screening can be proved by means 



§ . of helioseismology. 



I. INTRODUCTION 



o 

■ The solar neutrino problem is so important that any aspect of solar, plasma and nuclear 
physics pertinent to it has to be deeply investigated before definitive conclusions can be 
drawn. In this respect, screening of the charges of the reacting nuclei due to free charges in 

^ ' the solar plasma is of some interest. 

■ The study of screened nuclear reaction rates was started with the pioneering work of 
Salpeter who discussed both the extreme cases of "weak" and "strong" screening, pro- 
viding suitable expressions for the screening factors 

fij ij ^plasma / {'^'^'}ij,bare • (l) 

The solar core is not far from the weak screening case, however it does not satisfy the usual 
conditions under which the weak screening approximation holds. This is the reason why the 
problem has been investigated by several authors, see e.g. 

Gruzinov and Bahcall calculated the electron density in the vicinity of fusing nuclei 
using the partial differential equation for the density matrix that is derived in quantum sta- 
tistical mechanics. Their numerical result agrees, within small uncertainties, with Salpeter's 
weak screening formula. Furthermore, Bahcall et al. recently provided several arguments 
that demonstrate the validity of the Salpeter formula near the solar center with insignificant 
errors. 
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The conclusions of Gruzinov and Bahcall M are not unanimously accepted. According to 



Shaviv and Shaviv ||l0[, the weak screening formula does not hold in the Sun. Some nuclear 
reactions are enhanched by the surrounding plasma whereas some others are suppressed. 
According to Tsytovitch and Bornatici ITTIJI^ a kinetic description of collective plasma 
effects results in a decrease of the thermonuclear reaction rates in contrast to the Salpeter's 
enhancement. 

We observe that solar models are built by using stellar evolutionary codes which include 
specific expressions for the the nuclear reaction rates. If one uses different formulas for the 
screening factors fij one obtains different solar models. On the other hand, helioseismol- 
ogy provides precise information on the sound speed profile and on the properties of the 
convective envelope, see e.g. |]l3l, which have to be reproduced by the correct solar model. 
The main purpose of this paper is to test the screening models by means of helioseismology. 
We build several solar models corresponding to different screening factors and compare the 
results with helioseismic data. 

We also comment on the predicted neutrino fluxes. Many of the attempts to modify the 
screening factors have been produced as an effort to avoid or mitigate the so called solar 
neutrino puzzle, by reducing the predicted neutrino flux. We will show that a reduction 
of the screening factors does not generally imply a reduction of the neutrino flux. 



II. RESULTS OF SOLAR MODEL CALCULATIONS FOR DIFFERENT 

SCREENING PRESCRIPTIONS 



By using FRANEC [Ti], a stellar evolutionary code including diffusion of helium and 



heavy elements [0 , we constructed solar models based on four different assumptions: 



i) The weak screening approximation (WES). The screening factors fij are given by: 



In f^''^ = Z,Z,ey {an kT) (2) 

where Zi,Zj are the charges of the interacting nuclei, T is the temperature and is the 
Debye radius. As clear from equation above, the screening factors are always larger than 
unity, i.e. the plasma provides enhancement of the thermonuclear reaction rates. 

ii) The Mitler result ||^ (MIT), obtained with an analytical method which goes beyond 
the linearized approach and which correctly reproduces both the limits of weak and strong 
screening. Neglecting the small effects of a radial dependence in the effective potential, see 
0, the enhancement factors are given now by: 

In/f ^ = -linen,a%r [(0 + Q + if ' {Q + if - (0 + lf]/ikT) (3) 
where = SZi j/Airnea^^ and ng is the electron number density. 

iii) Neglect completely any screening effect (NOS), i.e. nuclear reactions occur with rates 
(cjf )bare- This casc is considered in connection with the suggestions that screening can be 
much smaller than Salpeter's estimate, see e.g. . 



iv) The Tsytovich model (TSY) [^,12], which provides a decrease of all the thermonuclear 
reaction rates with respect to the case of bare nuclei. Screening factors are taken from Table 
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lb of ref. and, for the ^Be electron capture, from Table 2 of the same reference (all 
factors are assumed constant along the solar profile). 

In Table we report the screening factors at the solar center for the various models. One 
sees that the weak screening approximation always yields the largest enhancement factors, as 
physically clear due to the fact that electrons and ions are assumed to be free and capable 
of following the reacting nuclei. We also remind that in this model the electron cloud is 
allowed to strongly condense around the nuclei. In the Mitler model, where electron density 
at the nuclear site is fixed at rif., the enhancement factor is smaller. By definition there is no 
enhancement in the NOS model, whereas in TSY model there is a decrease of the reaction 
rate, as already remarked. 

The main features of the solar models we obtained are presented in Table |l[ When 
moving from WES to solar models where nuclear reactions are less favoured one observes 
the following effects: 

i) The central temperature increases. In fact the hydrogen burning rate is fixed by the 
solar luminosity and a decrease of /n has to be compensated with a temperature increase; 

ii) The isothermal sound speed, u = P/p near the center increases. This is due to 
the increase of temperature whereas the "mean molecular weight" remains approximately 
constant; 

iii) The properties of the convective envelope are affected. In particular the border 
between the radiative and the convective region moves outwards and the photospheric helium 
abundance decreases. 



III. HELIOSEISMOLOGY AND ELECTRON SCREENING 



As well known several properties of the sun can be determined accurately by helioseis- 
mic data, see e.g. []T7|p!8| . The photospheric helium abundance Yp^ and the depth of the 



convective zone Rh are given by: 



Yph = 0.249(1 ± 1.4%) 



(4) 



Rb/Re = 0.711(1 ±0.2%) . 



(5) 

This error estimate 



The quoted errors are the so called "statistical" or " la" errors of |[I3|JT7 
was obtained by adding in quadrature each contribution to the uncertainty. Similar error 
estimates are given in [jl9|,^. A more conservative approach corresponds to add linearly 
all known individual uncertainties. This gives the so called "conservative" errors studied in 
||T3| , which are about a factor three larger than those in eqs. (§,^. 

Moreover, by inversion of helioseismic data one can determine the sound speed profile 
in the solar interior. This analysis can be performed either in terms of the isothermal 
squared sound speed, u = P/p, or in terms of the adiabatic squared sound speed = 
dP/dp\ad = iP/p, as the coefficient 7 = d log P/d log p\adiab is extremely well determined by 
the equation of state of the stellar plasma. The typical "Icr" error on u is about 1.3 °/oo in 
the intermediate solar region and increases up to 7 °/oo near the solar center, see [|T3|. A 
similar error estimate is obtained in 



18 . The "conservative" error estimate of ref 
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IS 



about a factor three larger. 
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Recent Standard Solar Models calculated by using the weak screening prescription are 
in agreement with helioseismic constraints on the properties of the convective envelope and 
on the sound speed profile, see e.g. the BP98 model of ref and BP2000 model of ref |^ 



As an example, we present in Fig. ^ the comparison between the prediction of BP2000 and 
helioseismic data for u = p/p- 

All this shows that the weak screening model is in agreeement with data and deviations 
from WES cannot be too large. From the comparison of different models, see Table || and 
Fig ^ one obtain the following results: 

i) The difference between the Tsytovitch model (TSY) and the weak screening model 
(WES) exceeds the "conservative" uncertainty on u in a significant portion of the solar 
profile. We remark that also the depth of the convective envelope is significantly altered. 
In other words the anti-screening predictions of ref. [|11],|12| can be excluded by means of 
helioseismology. 

ii) Also the difference between the no-screening model (NOS) and WES is significant for 
both u and Rb in comparison with helioseismic uncertainty. In other words the existence of 
a screening effect can be proved by means of helioseismology. 

iii) The Mitler model of screening (MIT) cannot be distinguished from the weak screening 
model within the present accuracy of helioseismology. 



IV. NEUTRINO FLUXES 

A reduction of the screening factors does not automatically mitigate the "solar neutrino 
problem". As an example the TSY model predicts a larger flux. Chlorine and Gallium 
signals than the WES model, see last column of Tab. 

As discussed extensively in ref 0, the behaviour of neutrino fluxes can be understood 
by considering that a decrease of the screening factors has the following effects on the solar 
structure: 

i) The hydrogen burning rate is flxed by the solar luminosity and a decrease of /n has to 
be compensated with a temperature increase, being approximatively @]: 

Tc oc /n'/' ; (6) 

ii) The rate of the ^He +^ He reaction, which is responsible for the PP-II chain and for 
Beryllium neutrino production, is changed. This results both from the increase in central 
temperature and from the variations in the screening factors f^^ and /ss, see As a 
consequence, one expects a variation in the beryllium neutrino flux given by : 



^ J 34 „-10/8 ,rj\ 



1/2 
33 



iii) The neutrino flux in addition depends on the ratio of the proton to electron capture 
rates on ^Be: 



'^BOC—-^-f\^ . (8) 

•^^7 /33 



,1/2 
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These scaling laws account for the numerical results of Table || and provide an expla- 
nation for the neutrino flux increase of the TSY model. As clear from eq.(|D effects on 
the proton and electron capture almost compensate {fij^^ / fri^^ = 2.213; f^^^ / fe7^^ = 
2.166), and the increase is essentially due to the /{^'^ term, which corresponds to the tem- 
perature effect. 



V. CONCLUSIONS 

We recall here the main points of our discussion: 
i)The ant i- screening predictions of ref pJ] , p!2| can be excluded by means of helioseismology. 



since both u and Rf, are significantly altered. 

ii) We find that the a no-screening solar model is not completely consistent with helioseismic 
data on u and Rb, in other words the existence of a screening effect can be proved by means 
of helioseismology. 
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TABLES 



TABLE L 


Screening factors 


in solar center, for weak screening 


5 (WES) 


|], Mitler model (MIT) 


IQ], no screenin 


g (NOS) and Tsytovitch model (TSY) 










WES 


MIT 


NOS 


TSY 






1.049 


1.045 


1 


0.949 


'^He +3 He 




1.213 


1.176 


1 


0.814 


^He +4 He 




1.213 


1.176 


1 


0.810 


'^Be + p 




1.213 


1.171 


1 


0.542 


TABLE IL 


Comparison among solar models with different 


screeninj 


g factors. We show the 


fractional differences, (model -WES )/WES, 


for the photospheric 


helium abundance (Yph), depth 


of the convective envelope (Rt), 


central temperature (Tc), isothermal sound speed squared at the 


solar center {uc 


), neutrino fluxes (<!>,) and predicted signals for Chlorine (CI) and the Gallium (Ga) 


experiments. All variations are in per cent. 










MIT 




NOS 


TSY 


Yph 


-0.076 




-0.86 




-1.4 


Rb 


+ 0.037 




+0M 




+0.59 


Tc 


+ 0.45 




+ 0.54: 




+1.4 


Uc 


+ 0.10 




+1.0 




+1.4 


%P 


-F0.033 




+0.45 




-0.35 




-0.19 




-2.4 




-5.9 


^SB 


-2.7 




-12. 




+11 


CI 


-2.5 




-11. 




+9.7 


Ga 


-0.76 




-2.9 




+2.3 
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FIGURES 



FIG. 1. Comparison between the BP2000 model |22| and helioseismic data for u = P/p. The 
"statistical" and "conservative" helioseismic uncertainties [13| correspond to the dark and light 
areas respectively. 



FIG. 2. Comparison of different screening models with the WES model for u = P/p, same no- 
tation as in Table I. The "statistical" and "conservative" helioseismic uncertainties correspond 
to the dark and light areas. 
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